Compositional (version 5.5)

Multivariate t random values simulation on the simplex: Multivariate t random values simulation on the simplex

Description

Multivariate t random values simulation on the simplex.

Usage

rcompt(n, m, s, dof, type = "alr")

Arguments

n

The sample size, a numerical value.

m

The mean vector in \(R^d\).

s

The covariance matrix in \(R^d\).

dof

The degrees of freedom.

type

The alr (type = "alr") or the ilr (type = "ilr") is to be used for closing the Euclidean data onto the simplex.

Value

A matrix with the simulated data.

Details

The algorithm is straightforward, generate random values from a multivariate t distribution in \(R^d\) and brings the values to the simplex \(S^d\) using the inverse of a log-ratio transformation.

References

Aitchison J. (1986). The statistical analysis of compositional data. Chapman & Hall.

See Also

comp.den, rdiri, rcompnorm, rmvt

Examples

# NOT RUN {
x <- as.matrix(iris[, 1:2])
m <- Rfast::colmeans(x)
s <- var(x)
y <- rcompt(100, m, s, 10)
comp.den(y, dist = "t")
ternary(y)
# }